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RIGOROUS  DERIVATION  OF  GENERAL  DISPERSION 
RELATIONSHIP  FOR  GYROTRON 


I.  Introduction 

1-3 

As  a high  power  microwave  and  mm  wave  source  the  gyrotron 

is  gaining  much  interest  recently  for  its  applications  to  RF  heating 

of  fusion  plasma,  to  radar  transmitter  and  to  amplifier.  Gyrotron 

makes  use  of  the  interaction  between  transverse  electric  field  and 

motion  of  relativistic  beam  electrons  via  the  electron  cyclotron  maser 

instability.  This  instability  is  originated  by  the  asymetry  (bunching) 

on  the  electron  beam  distribution  over  gyration  orbit,  caused  by  either 

the  signal  electromagnetic  wave  (amplifier)  or  the  initial  perturbation 

in  electron  distribution  (oscillator),  which  enhances  the  magnitude 

of  EM  wave  at  the  expense  of  the  transverse  energy  of  beam  electrons,  and 

which  in  turn  further  fuels  the  bunching  mechanism. 

. 4-1 1 

Most  of  the  previous  theoretical  analyses  nave  been  carried 

out  within  the  framework  of  linearized  Vlasov-Maxwell  system  under 

the  tenuous  beam  assumptions,  and  some  of  the  results  have  been  shown 

12-14 

to  agree  well  with  those  of  the  experiments.  There  were,  however, 

some  limitations  of  previous  analyses,  in  that  they  are  limited  to 

the  azimuthally  symmetric  modes  and  to  the  specific  equilibrium  distribution 

function.  In  this  paper  we  extend  our  analysis  to  include  the  azimuthally 

asymmetric  mode  with  general  equilibrium  profiles.  We  further  extend 

the  geometry  of  the  system  from  simple  waveguide  to  coaxial  waveguide 

in  order  to  examine  the  contribution  of  the  center  rod  in  the  gyrotron. 

The  derivation  of  the  dispersion  relation  consists  of  two  major 

parts:  (1)  the  orbit  integral  to  calculate  the  perturbed  distribution 

function  and  (2)  the  dispersion  equation  itself  as  moment  integral 

of  the  perturbed  distribution  function.  In  the  first  part  (the  orbit 
Note:  Manuscript  submitted  May  16,  1979. 
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integral)  we  follow  closely  the  analysis  by  Chu,  which  makes  use 
of  the  Bessel  sumamtion  theorem,  and  in  the  second  part  (the  dispersion 
equation)  we  use  a new  approach  in  order  to  retain  arbitrariness  of 
the  beam  distribution  function. 

After  careful  examination  of  self-consistent  distribution  function, 
it  is  found  that  the  distribution  function  is  determined  entirely 
by  three  constants  of  electron  motion,  namely  y (total  relativistic 
factor),  p (axial  linear  momentum)  and  P (canonical  angular  momentum). 

Z w 

With  these  constants  as  independent  variables  for  the  distribution 
function,  we  can  obtain  the  dispersion  equation  for  general  equilibrium 
beam  profiles.  We  add  that,  during  the  course  of  the  orbit  integral, 
we  include  the  contribution  of  the  term  that  is  proportion  to  3fft/sP 
(see  eq  (33)),  which  is  previously  ignored  by  other  authors. 

In  this  paper  the  authors  attempted  to  present  as  much  mathematical 
detail  as  possible  for  future  reference  at  the  risk  of  repetition 
of  some  previous  works.  We  emphasize  that  this  work  has  been  primarily 
the  ground  work  of  our  current  and  further  developments  on  this  subject, 
such  as  the  effect  due  to  spread  in  the  energy  and  momenta,  and  due 
to  the  presence  of  inner  conductor. 

In  Chapter  II,  the  basic  equations  and  assumptions  are  described. 
The  orbit  integral  for  the  perturbed  distribution  function  is  obtained 
in  Chapter  III.  In  Chapter  IV  the  dispersion  equation  for  arbitrary 


II.  Basic  Equations  and  Assumptions 


The  basic  equations  which  govern  the  combined  system  of  the  rela- 
tivistic electron  beam  with  electric  charge  (-e)  and  mass  m and  the 
electromagnetic  wave  are  as  follows: 
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where  the  quantities  B Oc,  t),  E(x,  t),  f(x,  £,  t),  n(x,  t)  and  J(x,  t) 
represent  the  magnetic  field,  the  electric  field,  the  distribution 
function  of  the  electron  beam,  the  beam  density,  and  the  beam  current 
density,  respectively.  Here  x and  £ denote  the  coordinates  in  real 
and  momentum  spaces,  respectively.  The  electromagnetic  wave  (E  and  B) 
is  determined  by  the  Maxwell  equations  (1)  to  (4),  the  electron  beam 
by  the  relativistic  Vlasov  equation  (5),  and  the  source  equations  (6) 
to  (7)  combine  these  two  systems.  The  Maxwell  equations  (1)  to  (4) 
can  be  combined  into  two  wave  equations. 
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Generally,  different  boundary  conditions  are  imposed  on  B and  E in 
equations  (8)  and  (9),  so  that  there  exist  two  separate  branches  of 
the  electromagnetic  wave  modes,  namely  transverse  electric  (TE)  and 
transverse  magnetic  (TM)  modes,  which  satisfy  equations  (8)  and  (9) 
respectively.  In  this  chapter,  we  will  examine  TE  modes  in  the  cylindrical 
waveguide  which  are  commonly  used  in  gyrotron  devices.  Hence,  the 
cylindrical  coordinate  system  (r,  9,  z)  is  employed  throughout  this 
paper. 

The  geometry  in  this  paper  is  assumed  to  be  coaxial  cylinder 
whose  conductors  are  located  at  r = 5R  and  r = R (Fig.  1).  The 
parameter  6 is  the  ratio  of  inner  conductor  radius  to  that  of  outer 
conductor,  and  note  that  when  6=0,  the  geometry  becomes  single  cylin- 
drical. 

Present  analysis  is  carred  out  in  the  framework  of  the  linearized 
Maxwell-Vlasov  system,  where  any  quantity  4»(x»  t)  is  expanded  into 

Mx,  t)  = ^(r)  + exp[i(kz  + &0  - i»t)]  . (10) 

Here  'J'q  is  the  equilibrium  quantity,  whose  symmetry  properties 

3 3 1 

(viz.  yy  = 0 = y^)  are  explicitly  assumed,  and  is  the  small  perturbed 

quantity  which  is  already  Fourier  decomposed  into  the  integer  azimuthal 

mode  number  SL , and  the  axial  wave  number  k.  With  eq  (10),  TE  mode 

equation  (8)  becomes 
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The  boundary  conditions  to  be  imposed  on  eq  (11)  are 
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The  perturbed  current  density  J is  given  by 

J*  = - e Jd^£  v f*  (14) 

The  perturbed  distribution  function  f*  is  given  by  orbital  equation 
using  characteristic  method 


3£' 


f1  = e j°  a(At)  exp[i(kAz+HA0  - toAt)]  (E*  + ^ v'  x B^  ) . (15) 

— CO 

where  primed  (')  coordinates  represent  the  equilibrium  location  of 

the  single  electron  at  time  t = t',  and  At  = t'-t,  Az  = z'-z,  A9  = 6'  - 0. 

It  is  evident  that  the  TE  differential  equation  system  (eq.  11,  14,  1: 
is  self-consistent  for  E*  and  B*.  A fully  self-consistent  problem 
is  difficult  to  solve,  and  some  assumptions  should  be  made  to  overcome 
this  difficulty.  First,  we  assume  that  the  electron  beam  is  very 
tenuous.  Specifically,  we  assume 

v/yQ  <«1  » 06' 
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where  the  Budker  parameter  V is  defined  by 


and  the  number  of  electrons  per  unit  axial  length  by 

N = 21T  Jrdr  / d3£  fQ  . (18) 

In  Eq.  (16),  y^  is  the  total  relativistic  mass  factor.  In  presently 

available  gyrotron  devices,  the  Budker  parameter  V is  in  the  order 
-3 

of  10  , so  that  assumption  (16)  is  valid.  Noting  that  the  beam  correction 

terra  in  the  right  hand  side  of  the  TE  dispersion  equation  (11)  is 
proportional  to  V,  under  the  tenuous  beam  approximation  (eq  (16)  ) we 
conclude  that  the  zeroth  order  solution  of  TE  waveguide  is  free  wave 
guide  mode.  We  will  use  this  mode  to  derive  the  perturbed  distribution 
function  f*  (eq.  15).  With  the  boundary  condition  (13),  the  free 
waveguide  mode  is  given  by 
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Here  unsigned  integer  n denotes  the  radial  mode  number,  and  and 
represent  first  and  second  kind  Bessel  functions  of  order  A.  The 
amplitude  of  B3  is  already  normalized  to  1 • If  we  combine  eqs  (11), 
(14),  and  (15)  with  those  in  (19),  and  by  making  use  of  orthogonality 
of  Bessel  function  (see  Appendix  A),  the  dispersion  equation  for  TE,^ 
mode  results.  Namely, 
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where  the  perturbed  distribution  function  f is  given  by  the  orbit 
integral  (eq.  15),  and 

<«L  - Z\a>  ' <‘*L  - *2)  /2aln  ■ (2U) 

Next,  we  assume  that,  at  equilibrium, 


E°  = 0 


£ - (0,  0,  Bq)  , 


(22) 


where  B^  is  constant  axial  magnetic  field.  These  assumptions  are 
consistent  with  the  tenuous  beam  assumption  (eq.  16)  in  that  we  neglect 
the  self-fields.  Therefore,  the  equilibrium  particle  orbit  is  composition 
of  circular  gyration  motion  in  transverse  direction  with  constant 
longitudinal  speed.  The  orbit  of  electron  is  discussed  in  detail 
in  Appendix  B. 

Another  assumption  we  make  is  the  symmetry  of  the  problem  at  equilibrium 
in  0 - and  z-  direction,  which  is  previously  made  use  of  in  eq.  (10). 

This  symmetry  enables  us  to  use  three  constants  of  motion  in  describing 
the  equilibrium  distribution  function  f^. 
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We  also  assume  that  the  growth  rate  UK  («*  Imw)  is  small; 
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where  the  classical  cyclotron  frequency  is  given  by 
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This  assumption  leads  to  single  mode  calculation  and  non-linear  inter- 
actions between  the  different  modes  can  be  neglected. 

Under  these  assumptions,  we  will  find  the  dispersion  relationship 
by  eq  (21).  This  requires  finding  the  perturbed  distribution  function 
f*  in  eq  (15)  via  the  orbital  integral  in  the  first  place. 
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III.  Orbit  Integral 
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and  guiding  center  radius  are  defined  in  order  (see  Appendix  B6-B8) 
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Now  let  us  make  a canonical  transformation  in  the  momentum  space 
£ to  the  variables  Y,  and  Pg  in  eq  (26).  Then  the  integral  I has 
the  form  below: 
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With  the  substitution  of  the  perturbed  TE  fields  in  eq  (19), 
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where 

Ij  = exp(-iC  ‘ ) Z^^^r’ ) - exp(+iC ' ) Zj^fa^r ' )]  » 

Z2  ~ r>  Zl  (ainr,)  * and 

I3  = ir'  sin  Z^a^r')  ■ i rQ  sin  X'  Z^^r’)  (35) 

Here  ZJ  (x)  = ^ Z^(x),  and  a',  X'  and  £'  are  angles  of  the  triangle 

OGA'  in  Fig.  3. 

Two  things  deserve  to  be  mentioned  in  the  orbit  integral  eq  (24) 
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with  its  integrand  I in  eq  (33).  First,  we  note  that  -g~,  -g^-  , -jp- 
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are  constant  along  the  equilibrium  orbit  path  and  can  be  factored 

out  of  the  integral.  Due  to  this  reason  we  introduced  the  set 

(y,  Pz>  P0^  aB  our  independent  variables.  Secondly,  the  integrands 

Ij,  l2»  and  1^  are  functions  of  r' , which  is  in  turn  a function  of 

t.  The  orbit  equation  r'  = r'(At)  is  not  only  complicated,  but  makes 

it  extremely  difficult  to  perform  the  integral.  We  note  that  the 

orbit  is  best  described  in  angle  variations.  They  are  from  Fig.  3, 

0) 
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where  X'“X  = ^4,»  (38) 

a - o’  = A©  • (39) 

It  is,  therefore,  desirable  to  convert  r'  (At)  into  r'  (A cp ) through 
angle  variations.  This  can  be  done  with  summation  theorem  of  Bessel 
function  (Appendix  A4).  Use  is  made  of  results  of  the  theorem,  i.e. 
eqs  (A-9),  (A-7)  and  (A-10)  to  obtain 
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end  a£n = ~j£~  was  previously  defined.  Here  upper  (lower)  expressions 

are  for  those  electrons  with  the  guiding  center  radius  to  greater 

(less)  than  the  Lannor  radius  r , that  is,  with  the  canonical  angular 

L 

momentum  P less  (greater)  than  zero.  Note  that  when  6=0  (viz. 
simple  waveguide  geometry)  the  combination  Bessel  function  Z is  replaced 
by  the  first  order  Bessel  function  J,  and  the  expressions  remain  the 
same  regardless  of  Pq.  Now  that  the  At  dependence  in  the  orbit  integral 
appears  only  through  exponential  function  with  angle  variables  X' , 
a',  0',  we  can  easily  perform  the  integral.  The  perturbed  distribution 
function  f1  is  proportional  to  the  angle  integral.  That  is, 
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represents  the  Doppler  shifted  cyclotron  frequency  of  harmonic  number  6. 
With  eq  (44)  the  full  expression  for  f*  is  given  by 
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I 

l 

(48) 

V 

1 

— 9 0)  — 
me  ^ | 

wvw 

z.-.<xo)j:<xl> 

(u>-n8)  ’ 

A 

x. 

8 * w ■ " j 

1 

1 l 

(49) 

I = 
z 

wc  k ~y 

i 

: (.<*-&  1 * 

1 8 

1 1 
j 

J.-t<,0)Zi<XL) 

1 

| ZU<X0)J.<XL>  ) 

<0 i j;-*(vz.<v  * 


, Z o (x  )J' (x  ) 

J s— X.  0 s L | 

+ XL  ) [ 

r - uc  n 

! + Y 

(50) 

1 J n (x«)Z' (x  ) ' 

L (o)-n  ) J 

Again  the  upper  (lower)  expressions  are  for  Pq  < 0 (P^  > 0),  and  x^, 

Xj  are  defined  by  equation  (43).  The  angles  Ot  and  y are  shown  in 
Fig.  4. 

It  should  be  noted  here  that  another  wave  mode  number  s is  introduced 

to  overcome  mathematical  difficulty  which  turns  out  to  possess  physical 

meaning  of  magnetic  harmonic  number  (see  eq.  45).  Also  note  that 

the  response  of  the  perturbed  distribution  function  gives  rise  to 

resonance  of  u)  at  SI  with  the  dependence  of  (u >-fl  ) The  discreteness 
s 8 

of  with  respect  to  s makes  it  possible  to  treat  the  problems  of 
different  mode  numbers  individually  (see  eq.  46).  We  will,  therefore, 
designate  the  transverse  electric  (TE)  mode  with  TEg^  notation  where 
l denotes  the  azimuthal  mode  number,  n the  radial  mode  number,  and 
8 the  magnetic  harmonic  number  of  cyclotron  frequency. 
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' 


..  — • 


m 


It  is  interesting  to  observe  the  responses  of  fl  to  the  independent 


momentum  variables  Y,  p and  P . While  I and  I (eqs.  48  and 

Y z 


49) 


v-1 


are  proportional  to  (u)-Jl  ) , I_  has  a constant  term  with  respect 

to  (oj-08).  Furthermore,  for  the  azimuthally  symmetric  (£,=0)  mode, 

_1  _ *■ 

the  resonance  terra  (u)-fl  ) completely  disappears  in  I_.  We  conclude, 

8 0 

therefore,  the  contribution  of  Sf^/Sp^  to  mode  will  remain  significant 

only  for  the  azimuthally  asymmetric  mode  (ij*  0). 


IV.  Dispersion  Equation  for  TEa  Mode 

a/FIS 

The  dispersion  equation  for  TE^ns  mode  is  obtained  by  substituting 
the  perturbed  function  f*  with  f1  into  the  original  TE  equation  (21). 

X/  II 

Namely, 

m2  J2  2 2ire  „ 1 

-2  Hn  c Hn  N£n 

x /rdr  Jdr£  fM  Z^jCo^r)  Civ.+v0)  + ZUl(<V)(iVv0>]-  (51> 

From  Appendix  B [eqs  (B-9),  (B-15),  and  (B-7)’,  we  have 

i vr  + ve  = exp[+  i(a+x)l  , (52) 


2 

d3p  - — — -A—  dy  dp  dpn  . (53) 

“c  r0  L 81nX  z 0 

Here  the  angles  a,  x>  and  £ of  the  triangle  OGA  are  as  shown  in  Fig.  4. 
With  eqs  (52)  and  (53)  substituted  into  (47),  the  dispersion  eq.  (51) 
can  be  rewritten  as  follows: 

4 - K K 

c An  Jin 


r 1 + _ A 3f 

x r0  sinx  ; ' t 3y  * 3PZ  8 9P0 


where 

♦ — = Z^+1(ajnr)  exp  [ i(e+l)x  + i(8*l)ot]  . (5: 

* 

and  I 'a  are  given  in  eqs.  (48)  to  (50).  We  note  here  that  I's  and  f^ 
is  independent  of  r (see  discussion  in  Appendix  C),  so  that  the  terms 
in  f ] bracket  can  be  factored  out  of  r-integration.  Thus, 


t 


i 


to  ,2  2 

2 K Hn 

2/  , V & A 3f  - A 3 f A 3 - q 


ainNitn 


rorLSinX 


(4>  - <T) 


when  (=  aj,nrL^  Previ°usiy  defined.  The  r-integration  in  eq.  (56) 
is  difficult  to  perform  through  complicated  dependence  of  angles  X 
and  & on  r (see  Appendix  B),  and  again  we  invoke  Bessel  summation 
theorem  (Appendix  A)  to  overcome  this  difficulty.  That  is,  we  expand 
Zfol  (°tj  r)  in  (eq.  55)  using  running  index  s'  on  triangle  OGA 
in  Fig.  4.  By  eq.  (A-6), 


vi‘vr>  ■ (-1)  U 


Zs'-£(X0)  Js'+1(XL)  \ 
Js'-Jl(x0)  Zs'+1(XL)  ^ 


x exp  [-i(s'+l)X  - i 01+  l)a]  , (5/ 

where  upper  (lower)  expression  is  for  P«  <0  (P0  >0),  and  x.  = a r . 

w u u x,  n o 

With  the  expansion  (57),  the  terra  ($+  - $ ) in  eq.  (56)  becomes 


$+  - 4>~  = 2(-lT  Z 


■ I 


Ze'-4<x0)  J.'(xL) 


Wx0>  z;>  (V 


x exp  [i(s-s' )x)  » (: 
The  upper  (lower)  expression  is  for  Pg  < 0 (Pg  > 0)  in  eqs.  (58)  and 
(59).  Using  eq.  (58)  in  eq.  (56)  leads  to 


i 


I 


4-*2-\n2--**r  to  KK  .'I-*!. 

c “in  Nin 


- 3fn  - - 3f 


J8,-£(X0)Z8,(XL) 


A VJL/\  A uln  ~ UJL  ^ 

[IYl?+I*ir+I9  ^ V 


where  $ , stands  for  the  integral  byr,  i.e., 


rorL  sinX 


exp  li(s-s’))(] 


Since  the  angle  X in  Fig.  4 is  function  of  r for  given  Y,  pz  and  PQ 

(i.e.  r and  r ),  we  transform  X = X(r)  in  eq.  (60)  via 
O L 

2 2 2 
rr>  +rT  "r 

cos  x - -§rr — ’ <61 

o L 

and  noting 


r rTsmX 

O L 


we  can  rewrite  4>  , as  integral  over  angle  x » 

S 8 


$ , = j*  dx  exp  [i(s-s')x  ] 


It  is  easy  to  show  that 

A 

* , - 2tt5  , , (64) 

88 1 88'  ’ 

where  6 , is  Kronecker  6,  which  eliminates  summation  over  s'  except 

for  8 1 *8 . Finally,  with  eq.  (64),  the  dispersion  equation  (59)  becomes 
as  follows: 

<4  - X2  - ^T-  • - M Jdp,  Si P. 

c R2  M.  * ® 

in 


X(TrTT  *T. 


♦ T.  ♦ if.  > 


Here,  from  eq  (21a), 

"*»  = “tX  ■ * X - j2)  Zfc2^Xtn7~^2x4n2-Z27Zl2^xJln7^  • 


When  P < 0 (or  rn  > r ),  we  use  the  following  T.'s  in  eq.  (65); 
0 U la  X 


Gd  2 

C X, 


T,E«!  ‘t'Vi'VW'  w 


> ’ 


T*  = k“c  -f  lz.-i(*o>W  > wn  • 


x«  = \vi<xo)j;(xl>  t <67) 

* \ z.-t(xo)j;(xL)  i 1 * 4^v- 1 1 

when  Pg  > 0 (i.e.  r^  < r^),  Zg_^(xg)  an§  Jg(xL)  are  replaced  by  js_£(xq) 

and  Zg(xL)  respectively  in  eqs  (66)  and  (67).  The  dispersion  relation 

eq.  (65)  can  be  further  reduced  by  integration  by  part.  However, 

caution  must  be  taken  in  the  process  especially  when  the  distribution 

function  f-  has  spread  in  y or  p . In  that  case  the  integral  over 
u z 

p , for  example,  has  a form  of 
z 


z kp  3p 

(JL>  ~ Z ~ 8—  °*Z 


and  the  denominator  ((D-fl  ) has  a simple  pole  at  p = — ? (<i)-s— ^),  which 

8 Z iC  s 

contributes  to  residue  at  the  same  value  p , thereby  resulting  in 

z 

Landau  type  damping.  The  same  is  true  for  yintegration.  We  have 
to  bear  in  mind  to  include  this  residue  contribution  whenever  need 


may  arise.  We  will  here  ignore  this  contribution,  by  noting  that 
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we  can  always  add  the  residue  contribution  to  the  results  of  integration 
by  part.  By  integration  by  part,  and  assuming  that  vanishes  at 

the  limiting  values  of  integration,  we  get 

2 o x*  2 o*2  2 

»L  . k2  . ta  . OjS-  Mip 


3T  3T  3T 

X [f  ) 1 

1 0 <5y  3p  3pT  J 


The  various  derivatives  in  [ ] bracket  can  be  computed  by  using  the 
table  I given  in  Appendix  B.  Then  we  obtain  for  eq(69) 


“ k2  . /d Y /dPz  f<»ef Q(y,  Pz,  P0) 

c2  R2  u^R 

r 2 (n  2-c2k2)  -« 

x - 6 ^ (x  ,x  ^ o 

x si  0 1 (w-n  )2 

s 

+ (u>-fl  ) { ~y  ‘WVV  ~ $2  JJsHs)l(xO’XL)  } 


[ + V (X0>XL)  " B Hs4(x0’XL)  j 

Here  is  defined  previously  in  eq  (66),  and 

"si'vV5  [ • 

s(s2-»  2) 

•WvV*  1 -r — z.-i(vw 

L 

x 

* T„  ZWV  W > • 


- ' - >•  ■ Wr - '**>», * WXVw W 


9 2 2 
2s  -x  x 


U #<Xn  ,X,)  5 J <XT)  j'(xT)  - \ <S-Jl)2i 

8 £,  0 L 8 L s L ' X-  2 f 


z;-*(v 


- XL  { Z8-/X0)}2  1 


We  define  the  following  quantities: 


kp  so) 

Si  = — - + — - 
s ya  y 


_ x!n  _ Xfn  _ Vl 

x0  ~ v \ = ~r  CL  ’ s = — 


rL  [»2c2(t2-d  - pz2i  11 , 

c 


r0  5 - Pz2  - 2mu)cP0]  ** 

c 


The  expressions  given  in  eq  (71)  hold  when  Pg  < 0.  When  Pg  >0, 


Z »(x.)  and  J (x_  ) are  replaced  by  J „(x  ) and  Z (x,  ) respective!-- 
s-Ji  0 s L s-x  o 8 L 


as  before. 


Once  the  distribution  function  f A Y,  p , P_)  is  given,  the  eigen- 

U Z w 


frequency  0)(  = (^  + iU). ) with  its  real  frequency  and  its  growth 


rate  |ux  | can  be  determined  from  the  dispersion  relation  eq  (70)  with 


terms  defined  in  eqs  (71)  and  (72).  It  should  be  emphasized  that 


the  dispersion  relationship  eq  (70)  is,  unlike  previous  works,  valid 


for  the  completely  general  distribution  function  of  the  beam  electrons 


as  long  as  the  cylindrical  symmetry  exists  and  the  self-field  can 


be  neglected.  The  left  hand  side  of  the  dispersion  relation  (70) 


is  the  free  waveguide  equation  with  cutoff  frequency  — — , while 

K 


the  right  hand  side  is  the  correction  to  the  waveguide  equation  due 


to  the  presence  of  the  beam.  Although  the  right  hand  side  is  small 


and  proportional  to  V (see  eqs.  (17)  and  (18)),  it  is  this  term  indeed 


J 
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» 


which  gives  rise  to  the  instability.  Careful  examination  shows  that 

\-2 


• - y , 

while  the  term  proportional  to  is  responsible  for  instability, 


the  terms  proportional  to  (w-^)'1  and  (uj-SM0  are  stabilizing  ones. 


More  detailed  analysis  of  dispersion  equations  will  be  carried  out 

15 


in  a separate  paper. 


V.  Zero-Temperature  Beam  Dispersion  Relationship 

In  this  chapter  we  will  derive  the  complete  dispersion  relationship 
for  a specific  distribution  function,  given  by 


u>  N 

Eo " ~;C2~2-  My-y)  ^Pz“Pz)  «(W  * 

4 it  me  y 


(73) 


where  y , p , and  are  constants.  This  electron  distribution  * ^c- 
’ rz  0 


- 2 


tion  describes  the  electrons  which  have  constant  energy  ymc  , constant 

guiding  center  radii  Tq»  and  constant  Larraor  radii  r . Here  the  constants 

r , and  r are  defined  as  follows: 

L U 


V ^ ' "2'2^2  * l>  - O 15  • 

7o  5 SST  1 "2c2  (?'u  ‘ *1  ‘ ’ 


(74) 


Since  there  is  no  spread  in  y,  p or  P , this  distribution  is  often 

Z U 

called  a Zero-temperature  distribution  function.  It  can  be  easily 
verified  that 

2Tr/rdr/d3£  fQ  = N , (75) 

where  N is  the  total  number  of  electrons  per  unit  axial  length.  After 
some  mathematical  manipulation  with  eq.  (73),  the  Zero-temperature 
beam  dispersion  relation  is  derived  below. 


id  , x 2 

— - k2  - 4=- 


2v  x 


dn 


YNdnR 


H „(x-, 
el  O' 

V 

2 I C 

o 

(uj-fi  ) ' Y 

Qsd 

s 

.i  <v  v 

- e2 

(n  2 - c2k2) 

8 

(u>-nx)2 


(76) 
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Here  barred  quantities  are  those  which  are  obtained  after  y = y, 


p * p , P “ P_  are  inserted  in  eqs  (71)  and  (72).  It  is  this  dis- 
persion  relation  we  will  examine  extensively  both  analytically  and 
numerically,  the  results  of  which  will  be  reported  in  separate  papers. 


VI. 


Conclusions 


A detailed  procedure  of  obtaining  a general  dispersion  equation 
for  the  electron  cyclotron  maser  instability  in  the  waveguide  has 
been  presented  with  and  without  a center  rod.  The  analysis  has  been 
carried  out  in  the  framework  of  linearized  Vlasov-Maxwell  system  under 
the  tenuous  beam  assumption  for  the  azimuthally  asymmetric  mode  and 
for  the  aribrary  equilibrium  distribution  function. 

The  orbit  integral  for  the  calculation  of  the  perturbed  distribution 
function  is  solved  with  the  aid  of  the  3essel  summation  theorem  without 
neglecting  any  term,  and  the  moment  integral  for  the  dispersion  equation 
for  an  arbitrary  equilibrium  beam  profile  is  obtained  with  the  help 
of  self-consistent  properties  of  the  distribution  function,  namely 
the  fact  that  the  distribution  function  is  entirely  described  by  three 
constants  of  motion.  As  an  example,  the  closed  form  of  dispersion 
function  for  cold  electron  beam  is  given  by  eq  (76). 

The  resulting  general  dispersion  equation  ((69)  or  (70))  is  useful 
for  further  development  of  the  theory,  such  as  the  effect  of  the  spread 
in  energy  and  momenta.  The  analysis  of  the  cold  beam  dispersion  relation 
is  given  in  accompanying  memorandum  paper.  ^ 
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APPENDIX  A 


Some  Properties  of  Bessel  Function 


In  this  appendix  some  properties  of  general  Bessel  function, 

which  are  used  in  the  text  are  summarized.  General  Bessel  function 

Z (x)  is  defined  by 
m 

Z (x)  = C J (x)  + C Y (x), 
m l m z m 

where  C,  and  C_  are  constants,  and  J and  Y are  the  first  and  the 
12  mm 

second  kind  Bessel  function  of  integer  order  m with  x>0 

1)  Recurrence  formulae 

Z . (x)  + Z (x)  = — Z (x) 

m-1  m+1  x m 


Z . (x)  - Z . (x)  = 2 Z'  (x) 
m-1  m+1  m 


(A-l) 


Z . (x)  Z (x)  = Z'  (x) 

m-1  x ra  m 


Z’  (x)  = - Z (x)  - Z . (x) 
m x m m+1 


2)  Negative  integer  order 


Z (x)  = (-1)“  Z (x) 
-m  m 


3)  Asymptotic  behavior  for  small  arguments  (m^O) 


J <x>-  i,  <S>" 
m m!  2 

JQ(x)-  1 - (|)2 


, (m  + 0) 


(A-2) 


(A-3) 


V>  * - ^ <;>"  • <■  * °> 

Vx)  - i <f> 

4)  Summation  theorem 

In  a triangle  OGA  as  shown  in  Fig.  4,  first  let 


a = < GOA,  £ = <0AG,  y = <0GA,  and  r = OA,  r^  = OG,  r^  £ GA.  Then  we 


5 ♦ a + x - * . 


(a-;) 


t 


Bessel  summation  theorem  states  that,  when  (i.e.  a >^), 

then  for  real  a > 0 , 

expiifcoOZ^iar)  = s?_a)  Z^+a(arQ)  Jfl(arL)  exp  (isX)  (A-5) 

with  integer  s.  If  r^  < (i.e.£  < -^) , we  have 

exp(iiC)  Z^(ar)  * JL®  ZZ+s^*rL^  Js^ar0^  exP^isX^  (A-5‘) 

From  here  on,  we  designate  unprimed  (primed)  equations  for  rn>r  (rft<r  ). 

U L U L 

We  can  rewrite  eq.  (A-5)  by  substituting  (-&)  for  £,  and  (-s)  for  8, 
and  by  making  use  of  eqs.  (A-2)  and  (A-4).  Namely, 

Z^(ar)  = (-1)^  E Zg  £(ar0)  «MarL)  exp[isy+Hot  1.  (A-6) 

Similarly  by  substituting  -%  for  i,  and  -a+v  for  s in  eq.  (A-5’), 

Z (ar)  = (-1)'  E J a (*rn)  z (ar, ) exp[isx+iZa  1 (A-6‘) 

After  we  differentiate  eqs  (A-6)  and  A-6')  with  respect  to  a,  we  get 
rZA(ar)  - (-1)*  Z tr^U^)  J/ar^  + ^^(ar,,)!^^)) 

x explisx  ♦ ila  1 , (A-7) 

rZjj,(ar)  = ^ | ^Vs-J^V  Z»(arL)  + rL  Js-ZUr0)  ZsUrL)1 

x exp[isx  + iia  ] • (A-7') 

By  substituting  Z+l  for  S-,  and  a+1  for  s into  eqs  (A-6)  and  (A-6'), 

we  obtain  with  eq  (A-4) 
e:,-p[+i-;]  Z (ar) 

- Wr0>  Js±,  (arL)  <A'8) 

■ (-I)*"  | exp[isX  + iZa  1 

Js-i(ar0)  Zs  ± l(arL)  (A"8,) 

Therefore,  from  eqs  (A-8)  and  (A-8')  with  eq  (A-l),  we  find 


*s[exp(*i£)  Z4-1(ar)  - exp(+i£)  Z^+1<ar)] 

exp  f isx  + iZa  1 


(-1)*  £ 

8 


z-i(*ro>  j;  (“i> 


J .(ar.)  Z'  (ar  ) 
s-Z  0 a L 


(A-9) 

(A-9') 


When  we  substitute  (s  + 1)  for  s in  eqs  (A-6)  and  (A-6'),  then  by 

making  use  of  eq  (A-l),  the  results  are 

i sinX  ZA(ar)  ■ (-l)*"  | explisX  + i*-a  1 

<a-Z) 
ar 

o 

(s-Z) 

«r„ 


. (ar.)  J' (ar. ) 

+ -2_ 

Z'  .(ar.)  J (ar  ) 

(A-10) 

8-X.  0 

8 L 

arL 

s-Z  0 s L 

.(ar.)  Z' (ar. ) 

♦ -i- 

J'  .(ar  ) Z (ar. ) 

| (A-10*) 

s-x.  0 

8 L 

arL 

s-Z  0 s L 

5)  Orthogonality 

Consider  Z (x)  defined  by 

5 V*>  - Tpr^,)  V«> 

such  that 

V5V  ■ 

where  X.  > 0 is  obtained  from 
Zn  — 

Zi  (XZn>  = ° ‘ 


(A-ll) 


(A-12) 


(A-13) 


Z-Cx)  thus  defined  is  in  the  form  of  eq  (A-l)  so  that  all  of  previous  relations 

X# 


are  applicable.  If  we  introduce  new  variable  r 3uch  that 
X. 


X = 


Zn 


r=“Znr  • 


then  it  is  to  show  the  orthogonality  of  Z^(.a^nr)  as  follows: 


(A-14) 


l(rir  Zt(al„r)  V0l»'r)  " V6  ' 

Here  the  normalization  factor  N^n  is  given  by 


NZnE 


2«  2 [ 

Zn 


2 -Z2)  Z2(X.  ) - («*X  * - z")  Z?(«X.) 


2„  2 . 2.  _2/ 


Zn 


Zn 


Z Zn 


(A-15) 


(A-16) 


Here  §m'  is  Kronecker  6-notation. 
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APPENDIX  B 


Electron  Orbit  Under  Constant  Magnetic  Field 


In  this  section,  we  will  examine  a single  electron  orbit  with 
charge  (-  e)under  the  constant  axial  magnetic  field  (eq  22).  The 
electric  potential  $ is  assumed  to  be  zero  and  the  vector  potential 
for  constant  magnetic  field  A = (OjA^Ccr),  0)  due  to  the  symmetry 
in  0 and  z.  Then  Lagrangian  L of  the  relativistic  particle  is  given 
by 

L = -me2  [1  - i2  (r  2 + r2  0 2 + z 2)]  **  - J r0  AQ  (0-1) 

c 


where 


'0 


h BQr 


(B-2) 


Since  the  Lagrangian  L is  cyclic  in  0 and  and  time  independent, 
there  exist  three  constants  of  motion  associated.  Namely, 

5 P0  * Y®r2  © " f r A0  = rp0  - mt^r2  = const. 


* = P 

as  z 


ymz  = pz  = const. 


„ J 2 2 * 2 4 

H = »c  p + m c 


= Y tac  = const . 


If  we  let 


p = tm2c2(y2-l)  - p 

■L  Z 


21*5 


rL5 


(B-3) 

(B-4) 

(B-5) 


(B-6) 

(B-7) 


mu) 


r 2 0 i 

rrt  = I r. J 

0 - L mu) 


2P«  h 

] 

c 


and  let  the  angle  between  and  e0  as  5 as  shown  in  Fig.  4. 


sin  £ e + cos  F $ 
r ^ 0 


(B-8) 


(B-9) 
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then 


2 2 2 
p = - ijjnu)  (r  -r.  ) = nJW  rrT  cos  K - r 
0 c 0 L c L c 


That  is 


2 2 r 
r + rn  - 2rr  cos  s 
U L 


(B-10) 


(B-ll) 


This  cosine  equation  describes  the  triangle  shown  in  Fig.  3 with  OA  = r, 

OG  = rn,  GA  = r , and  the  constants  r.  and  r correspond  to  the  guiding 
0 L 0 L 

center  and  the  Larmor  radii.  Some  useful  properties  result  from  the 
triangle  OGA. 


0 r L 

— ■ ■—  SS  ■ - - ■ 1=  ■ t ■ ■ 

sins  sinX  sm« 


(B-12) 


2 2 2 
r +rL  ~r0 

cos £ = — J— > sinC 

L 


,2  2W  2 2.  h 

[(r‘rl  )(r2  -r  )j 


2 2 2 2 2W  2 2.  *5 

r0  +rL  ~r  , sinX  - 1(t  Ut2  -r  )J 


2r0rL 


2c0rL 


2 2 2 
r +r0  "rL 


, sin** 


rl  “ ^r0_rL^  > r2  “ r0  + ri 


[(r2-ri2)(rL2-r2)^  ** 


Note  that  rn  (P  , Y,p  ),  r (Y,p  ).  From  eq.  (B-9),  we  find  that  if 

w y Z L Z 

Pa  > 0 (PQ<  0),  then  r < r (r  > rT ) . That  is,  the  orbit  of  electrons 
“ u U L U L 

with  Pg  > 0 (Pg  < 0)  encompasses  (does  not  encompass)  the  origin  of 


For  future  references,  the  following  table  for  the  derivatives 
is  useful  (see  Table  I). 
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It  is  our  intention  to  use  the  set  (Y,p  »P-)as  the  independent 

momentum  space  variables  p,  provided  that  the  Jacobian  | p ^ |/  0 

1 0 

We  now  proceed  to  complete  the  Jacobian.  The  Cartesian  momentum 

coordinates  p , p can  be  shown  as  (from  Fig.  4) 
x y 

p (sin  £ cos  0 - cos  K sin  0) 


(B-13) 


p^  = pj_  (sin  £ sin  0 + cos  £ cos  0). 


Then  the  Jacobian  is  found  by 
3(p  »P  ) 3p  3p 

I JL*.  X.  I = _* 

1 S(!y,Pq)  1 3y  3P 


3p  3p  2 

»•  xi  me 


3Y 


3P, 


u> 


r0rL 


sinx 


(B-14) 


ie  wi0 

3 

Therefore  the  volume  element  in  the  momentum  space  d p is  now  written  as 
2 


d3  * E£l  _I_  J_  dY  dp,  dP 

9 u r_rT  sinX  * e 

— C U L 

by  making  use  of  the  table  I and  the  relations  (B-12) 


(B-15) 
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APPENDIX  C 


Self-consistent  Electron  Distribution  Function 

In  this  appendix  we  will  describe  characteristics  of  the  self- 
consistent  distribution  function  of  beam  electrons,  f^  under  the  influent.? 
of  constant  magnetic  field  B^.  Using  the  constants  of  motion,  Y , 
p and  P as  the  independent  variables  in  momentum  space,  the  equilibrium 

Z U 

distribution  function  f^  is  generally  written  in  the  cylindrical  coordinate 


system  as  follows: 


where 


fo  = fo(r*  Y’  V V’ 


(C-l  ' 


r = x . e , 
— r 


= [1  + 1 *S 


2 2 
m c 


(c-: 


Pz  5 £‘ez  ’ 


P = (x.e  )(j>.e  ) - h mo)  (x.e  ) 

0 — r 9 c — r 

Here,  the  symmetry  of  equilibrum  in  0-  and  z-  coordinates  (viz.  ■fe  =0“4: 

*0 

is  explictly  assumed.  It  is  our  objective  to  show  that  ——  = 0,  from 
equilibrium  Vlasov  equation.  The  time  independent  Vlasov  equation  is  given 
by 

3f«  - 

0 . (C-3; 


3f0  e 3f0 

- * Tx 

3f  0 8fQ 

Now  we  will  obtain  the  expressions  of  tj-j  and  — when  f^  is  given 


by  eq  (C— 1 ) . 


Ill 


3'o  3r  3£0  3y  8f0  3",  3f0  3P.  Sf0 

"Ji'u  -s?  * w * ~k  sr  * ~e;  jf: 

— — — — Z — 0 


(04) 


3fom!L!!o  + n!fo+3^  !!o.+  !!e  ^ 

3£  3£  3r  + 3£  dy  * 3£  3p^  3£  3Pg 

It  can  be  easily  shown  from  eq  (02)  that 


3_r  _ ~ 3y  _ "rz 

3x  er*  3x  3x 


(05) 


- m 0 re 
c r 


(06) 


ll  = o , il.JL- 

3J>  ' * ^c2 


_l£  = g _5  = r g 

3£  V a£  e 


(07) 


g ^ 3 a 

Here,  the  relations,  90  ®r  = ®q  and  "39  ®0  = ” ®r’  are  made  u8e  °^* 

By  substituting  eqs  (06)  and  (07)  into  (04)  and  (05),  respectively,  we 


obtain 


3fft  3fn  3fn 


3fo  _ ^ afo  . 3fo  + rg  3fo 

32  y»2C2  3^  **  3PZ  6 3pe 


(08) 


(09) 


With  eqs  (08)  and  (09),  the  Vlasov  equation  (03)  becomes 


(OlO) 
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Since  V =0,  the  equilibrium  constraint  on  the  equilibrium  distribution 
function  f^  specified  with  eq  (C-l)  now  states  that 

3f 

■ ^ ^ ■ 0 (C-ll ) 

The  condition  (C-ll)  is  expected  from  simple  argument.  We  note  Vlasov 
equation  is  the  statement  of  Liouville  theorem, 


(C-12) 


If  the  functional  dependence  of  f^  is  assumed  as  eq  (C-l),  then 


Jf0  3£0  dr 


» w.  3 f dp  3 f dP 

12  E 0 2l  * 0 21  0 . o _e 

dt  3r  dt  9y  dt  + 9p  dt  3p  dt 

z 0 


(C-13) 


As  pointed  out  in  Appendix  B,  however,  the  variables  Y,  Pz  and  P^  are 

constants  of  motion,  so  that  = ^z  = ^^9  = 0 . Then  Liouville 

1 dt  dt 

theorem  eq  (C-12)  now  becomes 

T?  Tt  " 0 - <c- 

which  is  identical  with  previous  result  eq  (C-10). 


(C-14) 


We  see  the  significance  of  the  equilibrium  condition  0/3r  = 0 in 
3fQ  3fQ  3fQ 

that  this  makes  — g^  , -g^-  and  -g^-  to  be  independent  of  r,  thereby 
independent  of  gyration  angle  X (see  chaps.  Ill  and  IV)  so  that  these 
terms  can  be  factored  out  of  the  orbit  integral  (eq.  24)  and  of  r integra- 
tion in  the  dispersion  (eq.  51).  The  constants  of  motion  (Y,  p , Pn)  as 

z 0 

independent  variables  were  thus  adopted  in  the  description  of  equilibrium 
distribution  function  of  beam  electrons. 
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TABLE  I 


Derivatives  of  quantities  (Q's)  in  the  first 
column  with  respect  to  y (second  column),  j-<2 
(third  column)  and  PD  (fourth  column) . This 
table  was  used  in  Appendix  B. 


L„  l - Cross  section  of  coaxial  gyratron.  The  co™^tors 
re* located  at  r=R  and  r=6R,  and  the  electron  be,,  n- 

9ined  between  two  conductors.  The  constant  equilibrium 
agnetic  field  is  in  ^ direction. 


* © 


,<0  o . Fiectron  orbit.  The  transverse  equilibrium  orbit 


